Confront holographic QCD with Regge trajectories of vectors and axial-vectors 



Song He 1 , Mei Huang 1 * 2 , Qi-Shu Yan 3 ' 4 , and Yi Yang 5 ' 6 
1 Institute of High Energy Physics, Chinese Academy of Sciences, Beijing, China 
2 Theoretical Physics Center for Science Facilities, Chinese Academy of Sciences, Beijing, China 
3 Department of Physics, University of Toronto, Toronto, Canada 
4 Department of Physics, National Tsing Hua University, Hsinchu, Taiwan 
5 Department of Electrophysics, National Chiao-Tung University, Hsinchu, Taiwan and 
6 Physics Division, National Center for Theoretical Sciences, Hsinchu, Taiwan 

(Dated: December 6, 2008) 

We derive the general 5-dimension metric structure of the Dp — Dq system in type II superstring 
theory, and demonstrate the physical meaning of the parameters characterizing the 5-dimension 
metric structure of the holographic QCD model by relating them to the parameters describing 
Regge trajectories. By matching the spectra of vector mesons pi with deformed Dp — Dq soft-wall 
model, we find that the spectra of vector mesons pi can be described very well in the soft-wall 
D3 — Dq model, i.e, AdSs soft-wall model. We then investigate how well the AdSs soft-wall model 
can describe the Regge trajectory of axial- vector mesons oi. We find that the constant component 
of the 5-dimension mass square of axial- vector mesons plays an efficient role to realize the chiral 
symmetry breaking in the vacuum, and a small negative z 4 correction in the 5-dimension mass 
square is helpful to realize the chiral symmetry restoration in high excitation states. 

PACS numbers: 11.25.Tq, ll.10.Kk, 11.25.Wx, 12.38.Cy 



I. INTRODUCTION 



Quantum Chromodynamics (QCD) has been accepted as the basic theory of describing strong interaction for more 
than 30 years. However, it is still a challenge to solve QCD in non-perturbative region where gauge interaction is 
strong. In the early 1970's, string theory was proposed to describe strong interacting particles [l|. Recently, the 
discovery of the gravity/gauge duality Q has revived the hope to understand QCD in strongly coupled region using 
string theory. The gravity/gauge, or anti-de Sitter/conformal field theory (AdS/CFT) correspondence provides a 
revolutionary method to tackle the problem of strongly coupled gauge theories. For a review of AdS/CFT, see 
[||. The string description of realistic QCD has not been successfully formulated yet. Many efforts are invested 
in searching for such a realistic description by using the "top-down" approach, i.e. by deriving holographic QCD 
from string theory Q, as well as by using the "bottom-up" approach, i.e. by examining possible holographic QCD 
models from experimental data [fj, [fj, [D, l8|, [9J . 

It is an essential and crucial point for the realistic holographic QCD model to reproduce Regge behavior. Regge 
behavior is a well-known feature of QCD [10(, and it was the commanding evidence for suggesting the string-like 
structure of hadrons. A general empirical expression for Regge trajectories can be cast as 

M^ s = a n n + a s S + b, (1) 

where n and S are the quantum number of high radial and spin excitations, respectively. The slope a n and as have 
dimension GeV 2 , and describe the mass square increase rate in radial excitations and spin excitations, respectively. 
In principle, a n is not necessarily the same as as, though a n = as can be taken as a good approximation by fitting 
experimental data [llj. The parameter b is the ground state mass square, and it is channel-dependent. 

Currently, one of the most successful models derived from string theory is the Sakai-Sugimoto model [l^] , which 
can describe spontaneously chiral symmetry breaking naturally, but fails to generate the linear Regge behavior. On 
the other hand, in the "bottom-up" approach, many efforts have been paid to generate the linear Regge behavior 
for meson spectra [l3|, [3, HH and for baryon spectra [l6j]. In Ref. [Hj], Karch, Katz, Son and Stephanov (KKSS) 
found that a z 2 dilaton field correction in the AdS§ background leads to the linear behavior Af 2 cx n. 

In this paper, we carefully study realizing Regge trajectories of vector and axial-vector mesons in holographic 
QCD model. We take the data of the radial and spin excitations of p and a from PDG2007 [l7|, which are listed 
in Table|U To describe Regge trajectories for both (pi, p 3 ) and (a\, a 3 ), we use the general formula Eq. ([1]). From 
the experimental data, the parameters of Regge trajectories can be determined by using the standard x 2 fit. The 
parameters for (pi,/?3) mesons and their correlations read 

< = +0.91 ± 0.23 / 1 \ 

a p s = +1.08 + 0.39 -0.82 1 (2) 
6* = -1.09 ± 1.18 V+0.43 -0.81 \J 
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TABLE I: Vector and axial- vector meson spectra (in GeV). 

The parameters for (ai, a^) mesons and their correlations read 

< = +0.81 ± 0.22 / 1 \ 

a% = +0.88 ± 0.39 -0.82 1 (3) 

b a = +0.13+ 1.17 V+0.43 -0.81 l) 

From the Regge trajectories of p and a, one can read the information of chiral symmetry breaking in the vacuum 
and asymptotic chiral symmetry restoration in highly excited meson states: 

Chiral symmetry breaking: It is known that chiral symmetry is spontaneously broken in the vacuum, thus the 
observed chiral partners are not degenerate. From the Regge trajectories of the chiral partners p and a, the chiral 
symmetry breaking in the vacuum is reflected by the difference of the ground state square-masses b p and b a . The 
difference is as large as 1 GeV 2 , i.e. b a - b p ~ M£ - M 2 pi ~ 1 GeV 2 . 

Asymptotic chiral symmetry restoration in highly excited states: It is noticed that the chiral partner becomes 
more and more degenerate in high excitations [18|] • For example, for radial excitations of p\ and a\ , the mass square 
difference of chiral partners M 2 X — M 2 X = 0.9947 GeV 2 at n = 1, and this difference decreases to 0.3227 GeV 2 at 
n — 6; For the radial excitations of p% and 03, the mass square difference of chiral partners M 2 3 — M 2 3 = 0.6408 GeV 2 
at n = 1, and it decreases to 0.2736 GeV 2 at n — 3. 

We will only focus on the realization of Regge trajectories of p\ and a\ in this paper. The paper is organized as 
following: After the introduction, we derive the general 5-dimension metric structure of the Dp — Dq system in type 
II superstring theory in Sec. HU By matching the spectra of vector mesons p\ with deformed Dp — Dq soft-wall 
model, in Sec. IHII we determine the 5-dimension metric that can describe the spectra of vector mesons p\. We then 
investigate how to describe the spectra of axial-vector mesons a\ in the same background of p\ in Sec. IIV1 At the 
end we give summary in SecfVl 

II. 5-DIMENSION METRIC STRUCTURE OF D p - D q SOFT- WALL MODEL 
A. Metric structure of D v — D q system 

In order to investigate the possible dual string theory for describing Regge behavior, we introduce the following 
Dp — Dq branes system in type II superstring theory. In the Dp — Dq system, the N c background _Dp-brane describe 
the effects of pure QCD theory, while the Nf probe Dg-brane is to accommodate the fundamental flavors which has 
been introduced by Karch and Katz [lj|. Such a practice is well motivated from string theory side. For example, 
in Sakai-Sugimoto model, p = 4 and q = 8. Low energy hadronic excitations are fields on probe branes which in 
the background determined by the background branes. 

First, we consider N c background £>p-branes in type II superstring theory. The near horizon solution in 10- 
dimension space-time is [20| 

ds 2 = h~^r} IMV dx> J/ dx v + hh (du 2 + u 2 dQ,l_ p ) , (4) 
where //, v = 0, • •■ ,p, the warp factor h (u) — {R/uf p and R is a constant 

1 

7-JJ 

(5) 

The dilaton field in this background has the form of 



R = 



-fV(5- 
7T 



7-p 



g s NJ 7 s - 



e* =g s h{u) — 



(6) 
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The effective coupling of the Yang-Mills theory is 



9eff ~ 9sN c u p 3 , 



(7) 



which is u dependent. This u dependence corresponds to the RG flow in the Yang-Mills theory, i.e. the effective 
g e ff coupling constant depends on the energy scale u. In the case of D3-brane, g e ff ~ gsN c becomes a constant 
and the dual Yang-Mills theory is JV = 4 SYM theory which is a conformal field theory. The curvature of the 
background (|4]) is 



K 



1 



l s9eff 



(8) 



which reflects the string/gauge duality - the string on a background of curvature 1Z is dual to a gauge theory with 
the effective coupling g e ff- To make the perturbation valid in the string side, we require that the curvature is 
small 1Z <C 1, which means that the effective coupling in the dual gauge theory is large g e ff ^> In the case of 
D3-brane, the curvature 1Z becomes a constant, and the background ^ reduces to a constant curvature spacetime 
-AdS 5 xS 5 . 

The coordinates transformation (for the cases of p ^ 5) 



P \ p ~ a n szl - 
Rp~ 5 zp- 5 , 



brings the above solution ^ to the following Poincare form, 



ds 2 = e 2A ^ 



r] l j lU dx fl dx 1 ' + dz 1 



(p-5) 2 



l dfli 



(9) 



(10) 



Next, we consider Nf probe Dg-branes with q — 4 of their dimensions in the S q 4 part of S 8 p , with the other 
dimensions in z and directions as given in Table [III The induced q + I dimensions metric on the probe branes is 
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TABLE II: Spacetime embedding of D p — D q system. 



given as 



ds 2 = e 2A 



T] IMV dx tJ 'dx 1 ' + dz 2 H — ^dfl 2 _ 4 

Zn 



Where /i, v = 0, • • • ,3, and the metric function of the warp factor only includes the logarithmic term 

p-7 



A{z) = —ciq lnz, with qq 



and the dilaton field part takes the form of 



2 (p-5)' 



(11) 



(12) 



(p-3)(p-7) 



b — p z 



(13) 



It follows that 



3>(z) ~ dp biz, with do 



(p-3)(p-7) 
2(p-5) 



(14) 
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B. The deformed D p — D q soft- wall model 



In the above subsection, we have derived the general metric structure of the D p — D q system in Type II superstring 
theory, and we have noticed that the metric function A(z) only includes the logarithmic term, and in general there 
is another logarithmic contribution to the dilaton field. However, from the lessons of AdSs metric (D 3 system) and 
the Sakai-Sugimoto model(£>4 — D$ system), the D p — D q system cannot describe linear trajectories of mesons. It 
was shown in Ref . [l3j . in order to produce linear trajectories, there should be a z 2 term, but all z 2 asymptotics 
should be kept in the dilaton field &(z) and not in the warp factor A(z). Otherwise, the radial slope a n will be spin 
dependent. Therefore, to describe the real QCD, we propose a deformed D p — D q soft-wall model which is defined 
as 

A(z) = —ao lnz, $(z) = do lnz + c^z 2 . (15) 

By assuming that the gauge fields are independent of the internal space S q ~ 4 , after integrating out S q ~ A , up to 
the quadratic terms and following the same assumption as in fl3j . we can have the effective 5D action for higher 
spin mesons described by tensor fields as 

I= l -J d 5 x^ e-*w{A^ Ml ...M s A i V Ml '-- Ms 
+ m^ Ml ... Ms ^- M 4, (16) 

where 4>m 1 ---m s is the tensor field and Mj is the tensor index. The value of S is equal to the spin of the field. The 
parameter m 2 is the 5D mass square of the bulk fields, and g and $(z) are the induced q + 1 dimension metric and 
dilaton field as shown in Eq. (Ill|) and (|14|) . The action for pi, a\ and p3, 03 mesons is given by taking S = 1, and 
5 = 3 respectively. 

Following the standard procedure of dimensional reduction, we can decompose the bulk field into its 4d compo- 
nents and their fifth profiles as <fr(x; z)m 1 --M s — J2 n =o 4>m 1 ---Ms i x )' i l } n{ z )- The equation of motion (EOM) of the 
fifth profile wavefunctions ip(z) n for the general higher spin field can be derived as 

d 2 ^ n - d z B ■ d z ^ n + (Mis - m\e 2A ) ip n = , (17) 

where M n ^s is the mass of the 4-dimension field ^Mi - m s ( x ) anc ^ 

B = $-fc(2S , -l) J 4 = $ + c (2S'-l)lnz (18) 

is the linear combination of the metric background function and the dilaton field. The combination function B(z) 
approaches logarithmic asymptotic at UV brane, and goes to z 2 asymptotic at IR region. It is worthy of remark 
that the spin parameter S enters in the factor B and can affect the EOM and spectra. Two comments are in order: 
1) The EOM for the eigenspectrum and wavefunction is valid for all mesons. 2) The spin parameter S enters in the 
factor B and can affect the EOM and spectra. 

The parameter k is a parameter depending on the induced metric pip of the D q brane. After integrating out 
S q ~ 4 , k is determined as 

fc = _(p-3)( g -5 )+ 4_ 
7-p 

It is obviously that k depends on both p and q. For simplicity, we have defined 

(p- 3) fo-5) +4 

C0 = ka = - lP 2 y_ 5 / ■ (20) 

The parameters ao, do, the other two parameters k, cq and corresponding curvature 1Z for any Dp — Dq system 
are listed in Table Hill We notice that do = for D3 background branes, i.e. dilaton field is constant in AdSs space. 
However, the dilaton field in a general Dp — Dq system can have a lnz term contribution, e.g. in the DA — D8 
system (Sakai-Sugimoto model ,12]), do = —3/2. We also want to point out that for pure Dp — Dq system, the 
curvature is proportional to the inverse of the coupling strength g e ff- For D 3 background branes, the curvature is a 
constant. The curvature for D4 background branes is small at IR, and large at UV, its dual gauge theory is strongly 
coupled at IR and weakly coupled at UV, which is similar to QCD. However, the curvature for Dq background 
branes is large at IR, and small at UV, its dual gauge theory is weakly coulped at IR and strongly coupled at UV, 
which is opposite to QCD. 
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TABLE III: Theoretical results for the D p — D q system. 



III. DETERMINE HOLOGRAPHIC QCD MODEL FROM REGGE TRAJECTORIES OF p x 

In this section we firstly consider the spectra of vector mesons, and we will investigate axial-vectors mesons and 
the chiral symmetry breaking mechanism in more detail in Sec. IIV1 

In the dictionary of AdS/CFT, a /—form operator with conformal dimension A has 5-dimensional square-mass 
m§ = (A - /)(A + / - 4) in the bulk [Hj], and for vector mesons m\ = 0. The spectra of EOM of Eq. ([17} for 
vector mesons has an exact solution: 

Ml s = 4c 2 n + 4c 2 c S + 2c 2 (l - c + do) . (21) 

When Co = C2 = 1 and do = 0, this solution reduces to the results given in the Ref. [liij . 

This exact solution supports the parameterization on Regge trajectories and can tell how phenomenological 
parameters a n , as and b are directly related with the metric parameters cq, ci and d^: 

1) C2 is completely determined by the string tension in the radial direction a n , i.e. 

C2 = ±. (22) 

Andreev in Ref. [22j shows that there is an upper bound of c 2 . c? can be determined by the coefficient C2 [23| of 
the quadratic correction to the vector-vector current correlator [24| 

AV ^ = Co + 1(7 2 + £ -^2„«? 2n > , (23) 

n>2 

where C2 can be determined from e + e~ scattering data [2^|. According to [12], the relation between Ci in Eq. 
(|23|) and the parameter c 2 for Dilaton field is: C2 = —\C%. The experimental bound | Ci \ < 0.14 GeV 2 gives 
I C2 I < 0.21 GeV 2 . From the fitting result for the Regge trajectories of vector and axial-vector mesons Eqs. ([2]) and 
(|3|), we can read that c 2 is around 0.2 ~ 0.25, which is in agreement the experimental upper bounds. 

2) It is interesting to notice that Co reflects the difference of string tension in the radial direction and spin direction, 

co = 2* (24) 
a n 

From the string theory side, it is commonly believed that the dual string theory of describing QCD should be 
strongly curved at high energy scales and weakly curved at low energy scales [26j. It seems D4 background brane 
system is more like the dual string theory of QCD. However, by reading the result for the Regge trajectories of 
vector and axial- vector mesons Eqs. ^ and ([3]), we can see that as/a n is around 1, which indicates that the real 
holographic QCD model is not far away from AdS$ model. 

3) Furthermore, do can be solved out as 

dD= 26 +2 5_ 1> (25) 
an a n 

If we take the approximation of a n = as = 1, we have Co = l,c 2 = 1/4 for both vector and axial- vector mesons, 
while do is mainly determined by the ground state square-mass as d p J a — 2b p l a . 

From the above analysis, we can see that the realistic holographic QCD model should be close to models defined 
in Dp-branes background for p = 3, i.e., AdS$ background. By using Eq. (|2"Tj) to fit the spectra of vector mesons 
pi, we get the following metric parameters: 



1. 



% = 0.2GcV 2 , (26) 



c. 



G 



It is noticed that here, c 2 has dimension GeV 2 . In doing numerical simulation, we have introduced a dimensionless 
parameter u. The relation between u and z is defined as 

u = Ksci z . (27) 

In our fitting, we have fixed hsd = 0.2, uuv — 0.1 and ujjv = 3, respectively. We will use the same parameters for 
calculations of axial-vector mesons. 



IV. AXIAL VECTOR MESONS IN THE DEFORMED AdS 5 MODEL 

We have determined the metric structure Eq. (|26p of describing the vector mesons, and we suppose the spectra of 
axial-vector mesons ai can be described in the same holographic QCD model. As mentioned in the introduction, in 
order to describe the Regge trajectories of axial-vector mesons a, we have to describe the chiral symmetry breaking 
in the vacuum and the asymptotic chiral symmetry restoration in highly excited meson states. 



A. Chiral symmetry breaking 

In order to produce the splitting of axial-vector meson spectra from vector meson spectra, it is essential to know 
the chiral symmetry breaking mechanism, in the following we discuss two different ways of chiral symmetry breaking. 



Constant m\^ a in the Higgsless model 



If we start from the Higgsless model Eq. (fl6|). the only difference between the EOM for p and a is the 5D bulk 
mass m|. We have taken m| = for vector mesons, we can assume that m\ ^ for axial- vector mesons due to 
chiral symmetry breaking. 

For the general A(z) and $(2) parameterized as 

A(z) = -colnz, (28) 
*(z) = c 2 z\ (29) 

the potential derived from B s (z) — <&(z) — (2S — l)A(z) can be put as 

2 1 9 

777 — — 111 

V(z) = 2c 2 (m-l) + — T ^ + c 2 2 z 2 + ^ ; , (30) 

with 

m=[(25-l)co + l]/2. (31) 
When Co = 1, for any C2 and constant TO5, the EOM 

- + V(z)* = ml * . (32) 

can have analytic solution: 



m 2 n = c 2 (4n + 1) + c 2 (2S - 1) + c 2 ^[{2S - 1) + l] 2 + Am\ . (33) 

When 5 = 1, 



ml = 4c 2 n + 4c 2 + 2c 2 1^1 + m\ - 1 J . (34) 
From Eq. (j3~4"| . for m\ = and m\ a 7^ 0, we can get the spectra of axial- vector ai through shifting the spectra 



of pi by 2c 2 (y 1 + m| a — 1) upward. In Fig. [TJ we explicitly show our numerical results of the spectra of axial-vector 

mesons a\ by solving the EOM of a\ Eq. (j32|) for different m\. It is found that when m\ a = 0.5GeV 2 , the produced 
a\ spectra agrees well with the experimental data. 
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FIG. 1: The spectra of axial-vector mesons ai by solving the EOM of ai Eq. 132]) for different m\. 

2. Higgs Mechanism 

In Ref. [TH, it is suggested that the axial field picks up a z-dependent 5d mass via the Higgs mechanism from 
the background scalar X that encodes the chiral symmetry breaking. The axial vector meson mode equation reads: 

d z (e-^e A ^8 z a n ) + [m 2 n - g 2 e 2A ^ X(z) 2 }e-^e A ^a n (z) = , (35) 

with <7g = I2ir 2 /N c . The linearized equation of motion for the scalar field X reads: 

8 Z (e-^e 3A ^d z X(z)^ + ie-^e bA ^X{z) = . (36) 

It is quite complicated to solve the coupled Eqs. (|35j) and ([36]) . We only discuss two asymptotic solutions of scalar 
field at UV (z = 0) and IR, respectively. 

The asymptotic form at UV has the form of 

X (z) z =f I mgZ + Ivjz 3 , (37) 

here the coefficient m q is the UV (z — 0) boundary condition given by the quark mass matrix, while the coefficient 
£ is the chiral condensate, which can be determined dynamically by the boundary condition in the IR. We show 
the produced a\ spectra from Eq. ([37)) in Fig. [2] when N c ~ 3. Fig. [2ja) is for different values of m q with fixed 
S = (400MeV) 3 , and Fig. [2jb) is for different values of £ with fixed m q — 3MeV. It can be clearly seen that the 
asymptotic form at UV introduces nonlinearity of the n dependence of M 2 , the z 3 term contributes more on the 
non-linear behavior. 

For large z (in the IR) on the background <& = c 2 z 2 the equation for X becomes 

X" - 2c 2 zX' + %X = (z>l), (38) 
z z 

the scalar field has a solution goes to a constant in the IR, i.e., X(z) — > const as z — > oo. X(z) = const means 
the effective 5D mass m 2 a — g 2 X(z) 2 is a constant. Then the spectra of a\ will behave the same as that in Sec. 
IIVA1I 

B. Chiral Symmetry restoration 

Though the slopes in radial direction for vector mesons a p n and axial-vector mesons a° can be roughly taken as 
the same, i.e, a p n ~ a^, in order to have chiral symmetry restoration, we need aP n > o°. To accommodate such a 




case, we can take the ansatz that the effective 5D mass m\(z) has the following form 

?lf = 5c\z 2 + S -^ + Sml (39) 

Substituting this ansatz into Eq. (|30|) . we can find the spectra can be exactly solved out as 

ml = 4c' 2 n + 2c 2 S + 2c' 2 S' + 5ml + 2 (c 2 + c 2 ) , (40) 

where c' 2 = \/ c 2 + 5c\ and S' = V S 2 + 5S 2 . Then equivalently, the chiral symmetry restoration demands that c' 2 
should be smaller than c 2 , which infers that 5c\ should be negative. It is not clear at the moment, how to get the 
negative correction 5c 2 z 4 in 5D mass square of axial- vector meson from Higgs mechanism, which deserves further 
careful study in the future. 



V. SUMMARY 



In summary, we have derived the general 5-dimension metric structure of the Dp—Dq system in type II superstring 
theory. We have shown the dependence of the metric parameters on the Regge trajectories parameters: The 
quadratic term in the dilaton background field is solely determined by the slope in the radial direction; The warp 
factor is mainly determined by the difference of the slope in the spin direction and the radial direction; The 
logarithmic term in the dilaton background field contributes to the ground state square-masses. 

It is commonly believed that the dual string theory of describing QCD should be strongly curved at high energy 
scales and weakly curved at low energy scales, it seems that the DA background brane system is more like the dual 
string theory of QCD. However, the ratio of the slope in the spin direction of the vector mesons over its slope in 
the radial direction as/a n is around 1, which indicates that the real holographic QCD model is not far away from 
AdSs model. 

We have shown that the spectra of axial vector mesons can also be described in the AdS$ soft-wall model, and 
a constant 5D bulk mass for axial-vector meson plays efficient role to realize the chiral symmetry breaking in 
the vacuum, and a small negative z 4 correction in the 5D mass square is helpful to realize the chiral symmetry 
restoration in high excitation states. 

The information in this study is important for a realistic holographic QCD model and our future study on 
the interactions of mesons (branching ratios and decay widths of mesons, their interactions, etc). In our current 
approach, we chose the Higgsless model to describe the radial and higher spin excitations of vector and axial- 
vector mesons, the pseudoscalar n is the zero mode of axial-vector field and there is no its radial and higher spin 
excitations. We leave the study of the Regge trajectories of scalar and pseudo-scalar mesons by using the global 
symmetry breaking model in our future works. 
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